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ANALYTICAL  AND  EXPERIMENTAL  INVESTIGATION  OF  PERFORMANCE 
OF  SUPERSONIC  EJECTOR  NOZZLE* 


Shen  Hui-11  and  Chen  Zhong-qing 
(Northwestern  Polytecbnical  University) 

SUMMARY 

This  article  introduces  a  theoretical  method  of  calculation 
for  the  flow  field  and  performance  of  a  supersonic  ejector  nozzle. 

It  Includes  such  calculations  as  those  listed  here,  that  is,  the 
real  sonic  line  at  the  exit  of  the  primary  jet  nozzle,  the  non- 
viscid  flow  flow  field,  the  adjustments  for  the  Influence  of  vis¬ 
cosity,  pumping  performance  and  thrust  performance. 

In  order  to  make  the  results  of  calculations  agree  even  better 
with  the  data  obtained  from  experiments,  one  can  employ  the  real 
sonic  line  and  not  the  planar  sonic  line  as  the  initial  base  line 
for  calculations.  The  real  sonic  line  Is  obtained  from  the  points 
of  intersection  between  the  constant  flow  angle  lines  in  the  throat 
region  of  the  primary  Jet  nozzle  and  the  Mach  lines  which  occur  in 
the  lip  region  of  the  primary  jet  nozzle.  First,  one  needs  to  cal¬ 
culate  the  non-viscid  flow  flow  field  of  the  primary  Jet  nozzle  and 
then,  after  that  is  done,  make  adjustments  taking  into  consideration 
the  effects  of  viscosity.  The  method  which  is  employed  in  this  arti¬ 
cle  for  adjustment  to  correct  for  the  effects  of  viscosity  consists 
of  replacing  the  original  geometric  coordinates  of  the  ejector  shroud 
with  the  corrected  geometrical  coordinates  for  the  shroud.  The  for¬ 
mer  are  obtained  by  taking  the  displacement  thicknesses  for  the 
mixed  region  and  the  boundary  layers  and  superimposing  them  on  the 
original  geometrical  coordinates  for  the  shroud.  Calculations  of  the 
actual  flow  field  of  the  jet  nozzle  and  the  pumping  performance  as 
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well  as  other  related  quantities  are  done  on  the  basis  of  the  cor¬ 
rected  coordinates.  The  method  for  making  these  corrections  is 
relatively  simple  and  adequately  accurate. 

These  calculations  were  done  on  a  "320”  digital  computer  and 
model  testing  was  done  at  a  ground  facility.  The  theoretical  and 
actual  results  are  in  basic  agreement. 

ILLUSTRATIONS  OP  PLOW  FORMS  [1,2] 

l.Higb  harmonic  flows  (Figure  1).  When  the  amount  of  flow  involved 
in  harmonic  flows  is  small,  the  main  flow,  after  it  flows  out  of  the 
main  jet  nozzle,  continues  to  increase  its  acceleration  in  the  ex¬ 
pansion  section  of  the  "flow  body  wall  surface".  The  harmonic  flow 
involved  accelerates  in  the  passage  formed  by  the  surface  of  the 
shroud  and  the  "flow  body  wall  surface".  Depending  on  the  dimen¬ 
sions  of  the  ejector  nozzle  involved  and  the  conditions  of  the  flow 
being  considered,  there  are  three  types  of  possible  situations. 

These  are  the  maintaining  of  subsonic  speed  in  the  harmonic  flow, 
the  existence  of  sonic  speeds  at  the  end  surface  of  the  jet  nozzle 
within  the  harmonic  flow,  and  the  existence  of  supersonic  speeds  in 
the  harmonic  flow  on  the  end  surface  of  the  Jet  nozzle.  Because  of 
the  effects  of  the  viscosity  of  the  gases  involved,  there  is  formed 
in  the  interval  between  the  main  flow  and  the  harmonic  flow  a  non- 
isobaric  mixed  region  and  there  is  also  formed  a  boundary  layer  on 
the  wall  surface  of  the  shroud.  However,  in  the  case  in  which  one 
is  dealing  with  a  high  harmonic  flow,  the  influence  of  viscosity  is 
relatively  small. 

2.  Low  harmonic  flows  (Figure  2).  When  the  quantity  of  flow 
involved  in  a  situation  is  relatively  small  Cor  approaches  zero), 
then  the  main  flow  takes  the  form  of  a  free  Jet  and  expands  into  an 
isobaric  area.  This  jet  adheres  to  the  wall  surface  of  the  expan¬ 
sion  section  of  the  harmonic  flow  and,  going  through  a  violent  recom¬ 
pression,  forms  a  shock  wave.  In  the  upper  reaches  of  this  section. 
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there  is  formed  a  dead  zone  by  the  low  speed  isobaric  gas  body  which 
is  adjacent  to  the  boundaries  of  the  isobaric  area.  Along  the  bound¬ 
aries  of  the  main  flow,  due  to  the  effects  of  the  viscosity  of  the 
gases,  there  is  also  a  mixed  region  which  separates  the  main  flow 
from  the  dead  zone.  Inside  this  mixed  region,  there  is  an  exchange 
of  mass  and  energy  between  the  main  flow  and  the  harmonic  flows. 

Due  to  the  effects  of  viscosity,  the  main  flow  continuously  pulls 
out  a  certain  amount  of  the  harmonic  flow  from  the  dead  zone,  and 
this  section  of  the  harmonic  flow  which  is  pulled  away  is  then  re¬ 
placed  by  the  harmonic  flow  gas  supply  system. 


Figure  1.  An  illustration  of  the  flow  in  the  high  harmonic 
flow  within  an  ejector  nozzle. 

1 — harmonic;  2 — main 


Figure  2.  An  illustration  of  the  flow  in  the  low  harmonic 
flow  inside  an  ejector  nozzle 
1 — harmonic ;  2 — main 
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CALCULATION  METHODS 


1.  The  calculation  of  the  real  sonic  line  [3].  The  real 
sonic  line  is  determined  after  one  determines  the  angle  of  flow 
of  the  transonic  area  in  the  throat  section  of  the  main  Jet  nozzle 
and  the  Mach  lines  of  the  lip  of  the  main  Jet  nozzle  (Figure  3). 

If  we  make  the  assumptions  that  the  flow  involved  is  stable,  irro- 
tational,  isoentropic  and  two-dimensional,  then  this  type  of  flow 
can  be  described  by  using  a  flow  function  equation,  that  is 

( 1  -  (pVpW/pM*,.  +  t2( pVpW.V,/*1]*.,  ,  . 

+  [i-(p'/pW4>i^-o 

In  this  equation,  the  subscripts  are  the  partial  derivatives  chosen 
for  x  or  y.  The  input  definition  is  q  —  (**  +  0  —  arctg(i»/«,)  per¬ 

taining  to  the  speed  variable  and,  in  such  a  case.  Equation  (1) 
becomes 

a«pVp)*(a<r/d*) ]/a*  +  (pVp)(i/*)0  -  f/wvjdv  -  o  (2) 

If  we  again  use  the  definition  du  —  (p/p*)(</*/f)  and  its  transformed 
speed,  then  Equation  (2)  can  be  simplified  to  be 

&9/dtS  -M(M  wv/dv  -  o  (3 } 

In  this  equation  4(M)  -  (p*/p)jO  -  M1). 

If  we  make  the  "tangential  gas  assumption”  k(M)  =  1,  on  the 
basis  of  the  fact  that  the  speed  graph  relationships  can  be  simpli¬ 
fied  into  the  form  of  Cauchy-Riemann  equations,  it  is  then  possible 
to  solve  the  flow  field  involved  by  the  use  of  the  complex  variable 
method.  Compressible  and  non-compressible  flows  are  related  by  the 
equation  d*  —  (p/fXdq/q)  ,  and  after  we  do  an  Integration  of  this, 
we  obtain 

»  -  In [ 2(q/qmXMJy/t  -  A#L)/<  1  +  J  \  +  (q/qmyM'J{  \  -  M'J)))  (  4  ) 

In  thi3  equation,  qm  and  M^  represent  the  speed  parameter  and  M  num¬ 
ber  for  a  matching  configuration.  Because  of  tbe  fact  that  the 
Cauchy-Riemann  conditions  are  satisfied,  because  the  boundary  con¬ 
ditions  are  at  the  boundaries  being  studied,  it  is  possible  to  make 


use  of  a  set  of  original  and  combined  complex  fluctions  with  the 
conditions 

F(w  —  id )  »  In  [(cosh  coj  —  cosh  c0t)/(cosh  to*  —  !)] 

In  this  equation 

u>i  “  (10  —  id  —  ajj)-  •  */«;  M|  ”*  (u/  —  u>«)  •  */« 

If  we  carry  out  a  differentiation  of  Equation  (5),  then  we  obtain 
the  complex  speed 

dF /d(to  —  id)  —•  (*/<*)  rinh  «2/(co*h  »2  —  co*hto,) 

—  (*/«)  sinh  wi/(co*h  o»2  —  1)  ) 

Concerning  the  use  of  this  to  reach  solutions  in  situations  in 
which  there  is  a  finite  conical  jet  tube,  after  we  introduce  the 
complex  variable  * -*  j  +  *y  ,  then  it  becomes  possible  to  take  a 
graphic  solution  for  the  speed  and  transfer  it  to  the  physical  plane. 
When  this  is  done,  then 

dz  -  (l/c— ' '•ndF /d(u  -  id)  ]d(ca  -  id) 


In  this  equation,  the  last  quantity  with  the  horizontal  line  re¬ 
presents  the  complex  conjugate  value.  Concerning  the  integration 
of  Equation  C7),  it  is  possible  to  obtain  the  point  (x,y)  in  the 
speed  plane.  The  numberical  constant  in  the  equation  can  be  adjusted 
in  this  way  to  cause  the  position  of  the  jet  nozzle  lip  in  the 
speed  plane  a)  to  correspond  to  the  point  (.0,1)  on  the  physical 

plane.  What  is  particularly  useful  for  us  in  this  type  of  analysis 
is  the  distribution  of  the  flow  angles  of  the  throat  region.  The 
line  along  which  the  angles  of  flow  are  equal  is  called  the  "con¬ 
stant  flow  angle  line".  In  order  to  obtain  the  equation  for  the 
constant  flow  angle  line,  we  carry  out  a  numerical  value  integra¬ 
tion  for  values  of  w  from  w  to  u>.  in  the  range  in  which  the  flow 
angle  lies  between  0  and  a. 

One  makes  use  of  the  characteristic  curve  method  in  calculating 
the  Mach  line  for  the  lip  of  the  Jet  nozzle.  The  main  equations 
which  are  used  in  thi3  process  are  set  out  in  a  following  section. 
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In  the  lip  section  of  a  Jet  nozzle,  the  expansion  is  of  a  central 
wave  type  and  satisfies  the  Prandtl-Meyer  function.  Concerning 
the  fact  that  the  flow  in  the  lip  section  of  the  Jet  nozzle  turns 
through  a  certain  angle,  the  expansion  process  can  be  divided  into 
separate  finite  numerical  procedures.  The  first  point  on  the  sonic 
line  is  the  point  of  intersection  (Figure  3)  of  the  first  constant 
flow  angle  (the  angle  of  flow  is  a)  and  the  Mach  line  for  the  lip  of 
the  Jet  nozzle.  As  far  as  the  other  different  points  on  the  sonic 
line,  they  are  solved  for  in  an  analogous  way. 


Figure  3.  The  real  sonic  line 
1 — constant  flow  angle  line;  2 — sonic  line;  3 — Jet  boundary 

2.  The  calculation  of  the  flow  field  of  high  harmonic  non- 
viscous  flow  and  its  characteristics  (.Figure  4) 

Cl)  Assume  that  the  harmonic  flow  is  one-dimensional  and  iso- 
entropic  in  its  nature.  Assume  that  the  static  pressures  for  the 
two  sides  of  each  point  on  the  boundary  of  the  harmonic  and  main 
flows  are  equal  and  further  assume  that  the  main  flow  is  non-vis- 
cous  in  nature  and  is  constrained  by  the  characteristic  curve. 

C2)  Concerning  conditions  which  are  already  known,  the  initial 
conditions  concerning  the  main  flow  include  total  pressure,  total 
temperature  and  the  amount  of  flow.  These  conditions  also  include 
the  total  pressure  of  the  harmonic  flow  as  well  as  its  total 
temperature,  and  the  external  boundary  atmospheric  pressure.  The 
geometrical  dimensions  of  the  ejector  nozzle  are  also  included. 
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(3)  The  principal  formulae  used  are  [6] 


G  —  T*  (8) 

ft  “  arcsin(l/M)  (9) 

iJ  -  1  +  12/Of.  -  1)]  sinWO*  -  l)/(4  +  1)  •  9]  (10) 

dy/dx  “  tg(<£±^) ,  dy/dx  “  tg  qb  ffrj1  dt  — •  0 

(*P/pX<Sm*—  1/kM1)  +9<f>-  -( tia.4>/MR)»n  (11) 

(*p/p)(V  M*  —  1/*MJ)  -*<*>  —  -(sia  d>/MR)a§  (12) 


Figure  4.  A  computational  diagram  for  an  ejector  nozzle. 

1 — characteristic  curve;  2 — an  already-known  point;  3 — line  of  flow; 
4 — computation  point;  5 — -  axis  of  rotation 

(.4)  Computation  procedures.  First,  select  an  arbitrary  amount 
of  flow  for  the  harmonic  flow.  Then,  on  the  basis  of  the  condi¬ 
tion  of  mutual  equality  of  static  pressures  between  the  main  flow 
and  the  harmonic  flow,  make  a  precise  determination  of  the  expan¬ 
sion  of  the  main  flow  in  the  exhaust  of  the  main  jet  tube.  The 
precise  determination  of  the  points  on  the  boundary  lines  of  the 
two  flows  involved  is  done  by  the  use  of  the  successive  substitu¬ 
tion  method.  In  this  method,  as  far  as  the  static  pressure  of  the 
harmonic  flow  is  concerned  on  the  one  hand,  we  write  out  the  conti¬ 
nuity  equations  for  two  cross  sections  of  the  harmonic  flow  and,  on 
the  other  hand,  we  write  out  the  equation  for  the  characteristic 
curve  which  appears  in  the  vicinity  of  these  points.  These  calcu¬ 
lations  are  continued  straight  on  to  the  exhaust  of  the  jet  tube. 


If  one  Is  dealing  with  the  case  in  which  the  harmonic  flow  is  sub¬ 
sonic,  then  the  determining  condition  Is  that  the  static  pressure 
of  the  harmonic  flow  in  the  exhaust  of  the  jet  tube  is  equal  to 
the  surrounding  atmospheric  pressure.  When  the  harmonic  flow  in 
the  exhaust  plane  is  supersonic,  then  the  harmonic  flow  has  a 
throat  section  and  the  M  number  for  this  cross  section  is  1.  The 
determining  condition  at  this  time  is  the  identity  of  the  amount  of 
flow  which  is  selected  with  the  amount  of  flow  which  is  capable  of 
passing  through  the  throat  section  of  the  jet  tube. 

3.  The  adjustment  for  the  influence  of  viscosity  [3,4,7J. 

As  far  as  the  subject  of  the  influences  of  viscosity  is  concerned, 
reference  [4]  is  concerned  with  the  adjustment  of  the  area  of  the 
cross  seciton  of  the  throat  section  of  the  harmonic  flow  in  terms 
of  its  shape.  Reference  [3]  then  concerns  itself  with  adjustments 
for  the  Influence  of  the  boundary  lines  of  the  main  flow  and  the 
harmonic  flow  as  well  as  with  suc$  adjustments  for  the  influences 
of  the  boundary  layers  on  the  external  walls.  The  former  method  is 
simple;  however,  the  degree  of  accuracy  is  less  than  excellent.  The 
latter  method  of  adjustment  is  relatively  complicated.  What  this 
article  presents  as  the  best  solution  is  a  method  of  adjusting  the 
coordinates  of  the  shroud.  The  procedure  for  this  method  is  (1), 
on  the  basis  of  the  dimensions  of  the  original  ejector  nozzle  shroud, 
calculate  the  non-viscous  flow  field,  (2)  calculate  the  displace¬ 
ment  thickness  of  the  mixture  layer  between  the  main  flow  and  the 
harmonic  flow,  (3)  calculate  the  displacement  thickness  of  the 
boundary  layer  on  the  inner  walls  of  the  shroud,  (.4)  from  the  dimen¬ 
sions  of  the  harmonic  flow  passage  cross  section  which  was  calcu¬ 
lated  on  the  basis  of  the  non-viscous  flow  field,  subtract  the  dis¬ 
placement  thickness  of  the  boundary  layers  Involved,  then  add  the 
displacement  thickness  of  the  mixing  layer  and  get  the  dimensions 
of  the  harmonic  flow  passage  cross  section  after  it  has  been 
adjusted.  After  you  have  done  this,  adjust  the  coordinates  for  the 
wall  surfaces  of  the  shroud,  (_5)  on  the  basis  of  the  coordinates  of 
the  wall  surfaces  of  the  shroud  after  they  have  been  adjusted,  make 
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another  use  of  the  method  which  was  used  to  make  the  non-viscous 
calculations  and  use  it  this  time  to  calculate  the  entire  flow 
field.  The  details  of  the  calculations  are  as  follows: 

(a)  Concerning  the  calculation  of  the  displacement  thickness 
of  the  mixing  layer  of  the  harmonic  flow,  6£,  if  one  recognizes 
the  fact  that  the  mixing  layer  is  isobaric,  then  the  speed  diagram 
Fig.  5  for  the  interior  of  the  mixing  layer  can  be  obtained  from 
the  various  equations  below  [4] ,  that  is 

*-•/«..-  (1  +  qp»)/2  +  [(1  -  <p»)/2](l/«''J)([V',<tf)  (13) 

™  0  +  <T»)/2  +  1(1  ■“  <p»)/2  ]  erf  jj 

In  these  equations,  9>>  — «»/«,,  «»  is  the  speed  of  harmonic  flow,  u 

cL 

is  the  speed  of  the  main  flow  and  u - x  are  the  directional  com¬ 

ponents  of  velocity. 

erfn-(2/«w)  jV'Vf, 

7  — <r(y/x).  o  Is  the  equivalence 
parameter  and  It  is  obtained  from 
the  empirical  radius  formula,  y 
and  x  are  the  coordinates  in  the 
coordinate  system  within  the 
layer  of  mixing.  The  flow  speed 
at  the  point  at  which  y  =  0  is 
«  — (*.  +  *»)/2,  that  is,  <p  -  (1  +  <?>»)/2; 
n  “  the  equivalent  coordinate,  Flgure  ^  A  speed  chart  for 

that  Is  to  say  that  it  is  only  the  interior  of  the  mixing  layer 

necessary  for  the  value  of  y/x  to  be 

the  same  and  then  the  value  for  u/u  will  also  be  the  same.  That 

a 

is  the  same  thing  as  saying  that  the  speed  diagrams  for  the  various 
cross  sections  of  the  mixing  layer  are  all  similar. 

If  we  select  the  Prandtl  number  to  be  1,  and  ignore  dp/Qt , 
as  well  as  assuming  that  Cp  is  a  constant,-  then  it  is  possible,  from 
the  speed  diagrams,  to  directly  obtain  the  diagrams  of  temperature, 
that  is  to  say, 

9 


a  -  t*/t :  -  [1/(1  -  <p>)ii(T?/Tt)  -q>>  +  ( i  -  t;/t:)  -vi  (14) 

In  these  equations,  T*  is  the  total  temperature  inside  the  mixing 
layer.  T*  is  the  total  temperature  of  the  main  flow.  And  T£  is 
the  total  temperature  of  the  harmonic  flow. 

If  we  are  dealing  with  a  case  in  which  we  already  know  the 
speed,  temperature  and  static  pressure  of  the  various  points  inside 
the  mixing  layer  (the  static  pressure  is  equal  to  the  static  press¬ 
ures  of  the  main  and  harmonic  flows  on  the  boundaries  of  the  mixing 
layer),  then  it  is  possible  to  solve  for  the  displacement  thickness 
of  the  mixing  layer.  The  speed  diagrams  and  temperature  charts 
which  were  spoken  of  above  are  both  given  in  the  form  of  values  of 
x,  y  coordinates  in  the  internal  coordinate  system.  The  relation¬ 
ship  between  these  and  the  coordinates  of  the  reference  coordinate 
system  x,  y  is  Y  =  y  -  y  .  The  jet  boundary  lines  within  the  mixing 
layer  use  i  for  their  representation  (see  Figure  5).  As  far  as  the 
point  at  which  y  =  y  is  concerned,  it  is  best  to  take  the  mixing 
layer  and  divide  it  up  into  two  sections.  It  Is  best  that  the 
amount  of  flow  which  goes  through  the  range  form  to  nRb  be  equal 
to  the  amount  of  main  flow  which  enters  the  mixing  layer.  On  the 
basis  of  the  momentum  and  mass  of  the  mixing  layer  and  analysis  of 
these  two  quantities  in  light  of  the  principles  regarding  their  con¬ 
versation,  we  can  solve  for  that  Is 

/i (*ii)  “  Ui(mw)  —  /i(»j*»)]/(1  —  v>) 

ix(n)  -  10  -  c2)**»]/((r?/7?)  -  civil 

+  ('  Id  -  Cl)v/(A- clv^Un  , 

J,K»  (15) 

/i(n)  -  to  -  c'.)vI>i*.]/1(t:/t:)  -  civil 
+  [(i-ciV/Ol-cV)]'* 

J*lK> 

In  this  equation,  C  is  the  Clark  number.  Concerning  the  relation¬ 
ship  c1  —  M*/[2/(*  —  i)  +  MM  ,  Ym  Is  also  solved  for  on  the  basis  of 
the  principle  of  the  conservation  of  mass  and  momentum  and  is  given 
by  the  use  of  similar  coordinates,  that  is 

1 f •  “  Wm/ *  “  vk*  —  1 1/(  i  -  v ~  vJ 1  Cl 6 ) 
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According  to  the  concept  of  the  displacement  thickness  of 

boundary  layers,  one  gets  the  displacement  thickness  of  the  har- 

« 

monic  flow  of  the  mixed  region,  ^j^j^^gin  form 

-<TSl/z  -  -  UTf/TT)  -  Ci.rn/,(ijf)/[(l  -  -  n~  (-17) 

In  this  equation, 

<x  -  (<r/<r,)  •  <r,,  rflj  <r,  -  12  +  2 . 7 5 8 C.,/ v7  (T-  Ci, ; 

-  u.  +  (cWTr:)/(GtVT:)m i  +  (g*n/77)/(g,V rhu 

*/*x  -  (i  +  7'*)'o  -  -p»);  c.,  -  c\(\  -  'pi)/ic:a  -  <Vty  +  a  -  ci)i 

» 

(b)  The  calculation  of  the  displacement  thickness,  ^boundary 

of  the  boundary  layers  on  the  wall  surface  inside  the  shroud.  If 

« 

we  use  the  displacement  thickness  >  of  the  boundary  layers 

to  estimate  the  influence  of  the  boundary  layers  on  the  inner  walls 
of  the  shroud  on  the  amount  of  the  harmonic  flow,  then  the  proper¬ 
ties  of  the  boundary  layers  are  calculated  on  the  basis  of  an 
equivalent  length  and  this  equivalent  length  is  defined  as 

X  -  (Gy')-'  j'  (Gy‘)Jx  Cl8) 

In  this  equation  <7 —  (.«//(!+ 0.2M1)]4  ,  and  0  is  selected  as 

1.20  or  1.25,  depending  on  the  Reynolds  number.  The  Reynolds  number 
which  corresponds  to  the  equivalent  length  as  determined  by  the 
equation  above,  that  is,  R^  is  calculated  by  using  the  total  pressure 
and  total  temperature,  that  is  to  say 

K, -(V>*KTVr*]M» 0  +».2Ml)-ia  Mg) 


In  this  equation  r*  —  ').SlT*Rlft(T*)/p*] 


Concerning  the  use  of 


the  corresponding  Reynolds  number,  R_,  to  calculate  the  displacement 

*  x 

thickness  <Shniindary  and  the  momentum  thickness  0,  when  the  order  of 
Rx  is  10^,  then  we  use  the  equations  below  to  express  the  relation¬ 
ships  involved,  that  is 

8i  —  0.046jt(!  +  0.8A/{),-4,J?,-1'*  C2Q) 

0  -  0.036*(1  +  0 .1  Mi)-*  "/?,-•  »  C21) 


When  the  order  of  R  is  10,  then  use  the  equations  below  in 

A 

order  to  express  the  relationships  involved,  that  is 
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aS  —  0.028x(1.0  +  0.8  Miy"R.-t'a 
9  -  0.022x( t.O  +■  0.lMl)~,nR,~,w 


C22  ) 

(23) 

These  equations  above  are  accurate  for  k  =  1.4.  Due  to  the  fact 
that  we  recognize  that  the  static  pressure  inside  of  the  boundary 
layers  involved  is  a  constant,  when  we  calculate  the  pressure  of 
the  separation  surface  6  ,  it  is  possible  to  recognize  the  fact 
that  what  we  are  figuring  is  the  surface  pressure  and  this  pressure 
acts  as  the  basis  for  a  repeated  substitution.  Because  of  this,  as 
far  as  these  calculations  are  concerned,  it  is  possible  to  carry  them 
out,  in  their  cyclical  form,  using  a  standard  coordinate  form  which 
is  the  same  as  the  one  on  the  surface. 

4.  The  calculation  of  thrust  performance  [5].  Engine  thrust 
is  given  by  the  relationships  set  out  below,  that  is 

R  -  (C/g)V.  +■  p.A.  -  pHA.  -  (G/g)V,  -Re-  ( G/g)V .  (24  } 

In  this  equation  is  the  total  thrust,  that  is 

Fc  —  PhA , 

On  the  basis  of  the  theorems  relating  to  momentum,  the 
tion  above  can  be  changed  to  be 

Re  “  F#  -+■  F»  +  R,n  —  PhA, 

and 

F«  —  *(00-0  +  teuCyMDAj!  (jfcft  Mm  —  1)  (26) 

F*  —  p,  •  ^*,(1.0  +  '  (27) 

Rltt  —  (F*  •  2*y  'g  (<*>,«)  +  f,  •  *  •  ybp*MW  l  +  tg?(<t,H)<i*  (28) 

In  these  equations,  cf  is  the  coefficient  of  friction  of  the 
wall  surfaces  of  the  shroud.  cD  is  the  coefficient  of  flow  of  the 
main  Jet  tube.  cv  is  the  speed  coefficient  of  the  main  jet  tube. 

And,  y  Is  the  wall  surface  radius  of  the  shroud. 

When  one  is  considering  the  influence  of  viscosity,  the  coeffi¬ 
cients  for  the  speed  and  the  amount  of  flow  are  solved  for 


e  qua- 
(25) 
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respectively  by  the  use  of  the  equations  presented  below,  that  is 


e'o  “  2*y  (  cos<t>Jy  —  sin  2 *ydy 

f'«“  ^  j  2*y  (co t<fidy  —  sin  <^x)co4<^|^/|^ j  2«y  (  cot  <j>dy  —  sin  4>dx) 


The  integrations  in  these  equations  are  carried  out  along  the 
sonic  line.  If  one  is  considering  the  influence  of  viscosity,  then 
the  coefficients  of  speed  and  amount  of  flow  for  the  main  jet  tube 
have  formulae  which  are  respectively 

eu  —  (1.0  —  Re"*  *)  2*y  (co*$Jy  —■  sin  Inydy 

e,  —  (1.0  —  fcRe"*  *)  2*y  (cos ifidy  —  sin  d>dx) cos / 

2*y  (  cos  <txiy  —  sia</>4x)j 

In  these  equations,  k1  *  0.l85/cos  a;  k2  *  0.144 ;  4>  is  the 
included  angle  between  the  flow  speed  on  the  sonic  line  and  the  x 
axis,  y  is  the  vertical  coordinate  for  points  on  the  sonic  line, 
x  is  the  horizontal  coordinate  for  points  on  the  sonic  line.  And 
Rc  is  the  Reynolds  number 

A  COMPARISON  BETWEEN  THEORETICAL  CALCULATIONS  AND  EXPERIMENTAL 
RESULTS 


(.31) 

C32) 


Concerning  the  use  of  the  computational  methods  and  sequence 
which  are  presented  in  this  article,  we  used  them  to  calcuate  the 
ejector  nozzle  flow  field,  the  wall  surface  pressure  distribution 
and  the  pumping  performance  of  the  Jet  tube.  Moreover,  we  made  a 
comparison  of  these  calculated  results  with  the  results  from  exper¬ 
iments  involving  model  simulations  of  the  Jet  tubes  involved.  The 
results  of  this  comparison  are  as  shown  below. 

1.  The  flow  field  and  wall  surface  pressure  distribution  of 
a  jet  ejector  nozzle. 
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The  curve  network  for  the  performance  of  an  ejector  nozzle, 
as  it  was  calculated,  is  as  shown  in  Figure  6.  The  calculated 
values  for  the  pressure  distribution  on  the  wall  surfaces  of  the 
shroud  and  the  results  of  the  related  experiments  are  as  shown  in 
Figure  7. 

2.  The  pumping  performance  of  the  ejector  nozzle. 

As  Figure  8  shows,  curve  1  represents  the  results  of  calcula¬ 
tions  of  the  non-viscid  flow  and  curve  2  is  the  result  after  this 
original  curve  is  corrected  for  the  effects  of  viscosity.  Figure  8 
shows  that,  after  the  adjustments  for  viscosity  are  made,  the  theo¬ 
retical  and  experimental  results  are  basically  the  same.  Figure  8 
also  shows  that  when  one  is  working  with  a  situation  in  which  there 
is  a  high  harmonic  flow,  the  error  which  appears  between  the  non- 
viscous  flow  calculations  and  the  data  obtained  through  experimenta 
tion  does  not  exceed  10? . 


Figure  6.  A  graph  of  the  performance  curve  network  for  the  ejector 
nozzle  under  consideration. 

1— the  wall  surface  is  the  measured  value;  2 — harmonic;  3 — main 


Figure  7.  The  pressure  distribution  of  the  wall  surface 
of  the  shroud. 

1 — experimental  values;  2 — calculated  curve 


CONCLUSIONS 


1.  There  is  an  adequate  correlation  between  the  calculation 
methods  introduced  In  this  article  and  the  experimental  results 
obtained  for  the  same  values . 


2.  The  method  which  the  article  presents  for  viscosity  cor¬ 
rections  is  relatively  simple  and  is,  at  the  same  time,  adequately 


accurate . 


0.30 

*Z.0.2S 

STh 
0.20 


0.01  0.02  0.03  0.04  0.05  0.06  0.07  0.08 
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Figure  8  M  ? 


1 — experimental  values;  2 — calculated  values  for  non-viscous  flow; 
3 — calculated  values  for  viscous  flow 
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TABLE  OF  SYMBOLS 


^  flow  function 
p  density 
a  speed  of  sound 
u  axial  speed 

v  radial  speed  or  the  flow  speed 
in  the  exhaust  of  the  engine 

p  pressure 

p  Mach  angle 

flow  line  single  (the  included 
angle  between  the  flow  line 
involved  and  the  horizontal 
axis) 

9  the  included  angle  between  the 
characteristic  curve  and  the 
vertical  axis  or  the  boundary 
layer  momentum  thickness 

X  equivalent  length 

y  radial  distance  divided  by  the 
radius  of  the  jet  nozzle 

x  axial  distance  divided  by  the 
radius  of  the  Jet  nozzle 

oj  transform  speed 

a  hemipyramidal  angle  of  the 
exhaust  of  the  Jet  nozzle 

T  temperature 

G  amount  of  flow 

M  Mach  number 

X  coefficient  of  flow  speed 

k  specific  heat  ratio 

R  thrust  or  gas  constant  or 
radius 

A  area 

F  impuse  or  functional  relation 
n  coordinates  on  the  Mach  line 


C  coordinates  on  the  Mach  line 

3  constant  whose  selection  is 

dependent  on  the  Reynolds 
number 

Y*  a  viscosity  coefficient  repre^ 
sented  by  the  use  of  a 
blockage  coefficient 

a\  the  displacement  thickness  of 
boundary  layers 
(+  =  boundary) 

Subscripts 

a  pertaining  to  the  main  flow 
b  pertaining  to  a  harmonic  flow 
c  main  Jet  tube  exhaust 
e  cross  section  of  engine  exhaust 
H  boundary  conditions 
G  total 

0  engine  intake 
SH  shroud 

Superscripts 
*  blockage  coefficient 
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ANALYTICAL  AND  EXPERIMENTAL  INVESTIGATION  OF 
PERFORMANCE  OF  SUPERSONIC  EJECTOR  NOZZLE 

Shen  Hui-li,  Chen  Zhong-qing 
(Norik- 1? enter*  Polytethnieol  tTnirrrnty) 


Abstract 

This  paper  presents  an  analytical  method  for  calculating  the  flow  field  and  per¬ 
formance  of  supersonic  ejector  nozzle.  The  calculations  involve  the  renl  sonic  line  at 
the  exit  of  the  primary  nozzle,  the  inviscid  primary  flow  field,  the  correction  for  vis¬ 
cosity  effect  snd  the  pumping,  and  thrust  characteristics. 

In  order  to  bring  calculated  results  into  agreement  with  experimental  data,  the 
real  sonic  line,  instead  of  the  plane  sonic  line,  is  taken  as  the  initial  bnse  line  of  calcula¬ 
tion.  The  real  sonic  line  is  obtained  by  joining  the  points  of  intersection  of  constant 
flow  angle  lines  in  the  throat  region  with  Mach  lines  at  the  lip  of  the  primary  nozzle. 

First,  the  inviscid  primary  flow  field  of  the  nozzle  is  calculated  and  then  corrected 
to  account  for  the  viscosity  effect.  The  method  of  correction  for  viscosity  effect  prn- 
posed  in  this  paper  replaces  the  original  geometric  coordinates  «f  the  ejector  shroud 
with  corrected  geometric  coordinates,  which  are  obtained  by  superimposing  on  the 
original  geometric  coordinates  the  displacements  of  the  mixed  region  ami  the  boundary 
layer.  On  the  basis  of  th*  corrected  coordinates,  the  actual  primary  flow  field  and 
pumping  performance  of  the  nozzle  are  then  calculated.  The  proposed  method  proves 
to  be  quite  simple  snd  accurate. 

Calculations  were  performed  on  a  “JftJO”  digital  computer,  and  model  tests  on  a 
ground  test  facility.  The  analytical  and  experimental  results  arc  found  to  he  in  fairly 
satisfactory  agreement. 
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analytical  and  experimental  investigation 
op  supersonic  ejector  nozzle* 

Wu  Da,  Yang  Zai-ming  and  Zhang  Rung 
(.Beijing  Institute  of  Aeronautics  &  Astronautics) 

SUMMARY 

In  this  article,  we  introduce  some  test  results  concerning 
ejector  nozzles  and,  on  the  basis  of  these  results,  we  offer  a  sim¬ 
plified  model  of  flow  configuration  within  ejector  nozzles.  On 
the  basis  of  this  model,  we  based  approximation  formulae  with  which 
to  calculate  the  thrust  performance  within  these  nozzles.  Finally, 
we  carried  out  actual  test  calculations  and  took  the  results  of 
these  calculations  and  compared  them  to  actual  measured  values. 

The  results  of  this  process  showed  that  the  calculated  values  and 
the  measured  values  were  in  sufficient  agreement  with  each  other, 

1.  PREFACE 

Concerning  supersonic  ejector  nozzles,  because  of  the  fact  that 
their  structure  is  simple  and  compact,  in  situations  where  one  is 
dealing  with  large  drops  in  pressure,  these  nozzles  have  excellent 
characteristics.  Because  of  this,  they  are  widely  used  in  super¬ 
sonic  fighters,  bombers,  transport  planes  and  civilian  aircraft. 

Because  of  the  fact  that  the  requirements  of  various  types  of 
aircraft  are  not  the  same,  the  differences  in  structure  between 
different  types  of  ejector  nozzles  are  very  great.  If  we  speak  In 
terms  of  the  actual  types  of  hardware  involved,  a  small  number  of 
possible  structural  designs  Is  chosen,  after  careful  comparison  and 
selection  from  among  a  large  number  of  original  types  of  design 
possibilities,  after  these  possible  designs  have  gone  through  static 
state  testing  of  the  ejector  nozzles  Involved.  This  is  the  normal 
procedure.  If  one  operates  in  this  way,  one  finds  that,  in  such 


This  paper  was  read  at  the  2nd  Annual  All-China  Engineering  Thermo 
physics  Technology  Conference  in  Hangzhou  in  November  of  1978, 
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related  areas  as  design  plans,  the  machining  of  test  components  and 
the  organization  of  experiments,  the  amount  of  work  expended  is 
large,  the  development  periods  involved  are  long  and  the  procedure 
is  not  economical.  It  cannot  be  doubted  that,  if  one  were  able, 
in  the  initial  selection  of  a  design,  to  make  performance  esti¬ 
mates  of  various  different  types  of  design  plans,  in  order  to  weed 
out  those  designs  whose  performance  was  relatively  unsatisfactory, 
and  because  of  this  preselection  process,  only  have  to  organize 
evaluation  tests  for  a  small  number  of  designs,  that  the  process  of 
selecting  the  best  design  could  be  made  faster,  better  and  more 
economical . 

In  the  last  few  years  we  have,  on  a  floating  sleeve  type  test 
bed  for  the  evaluation  of  perpendicular  Intake  ejector  nozzles, 
made  tests  on  several  score  of  different  structural  models  for  ejec¬ 
tor  nozzles,  photographed  several  types  of  these,  as  well  as  the 
flow  spectra  for  these  nozzles  under  different  types  of  flow  condi¬ 
tions.  Moreover,  we  carried  out  analyses  of  these  experimental 
data  and  flow  spectra  and  came  up  with  types  of  models  for  the  inter¬ 
ior  of  ejector  nozzles  under  different  conditions  of  flow  configur¬ 
ation.  On  the  basis  of  these  models,  we  set  up  methods  for  the 
estimation  of  the  interior  performance  of  the  ejector  nozzles  being 
evaluated.  What  follows  is  a  discussion  of  the  main  experimental 
results  which  we  got  and  the  performance  analyses  we  did. 

2.  A  THEORETICAL  ANALYSIS  OF  THE  INTERIOR  CHARACTERISTICS  OF 
SUPERSONIC  EJECTOR  NOZZLES 

Cone)  Thrust  characteristics  for  the  interiors  of  ejector 
nozzles . 

On  the  basis  of  an  analysis  of  the  flow  spectra  of  the  ejector 
nozzle  exhaust  flow  of  ejectors  being  tested,  we  offer  a  type  of 
flow  model  for  main  Jet  nozzles  placed  in  configurations  of  incom¬ 
plete  (or  complete)  expansion,  and  we  make  the  following  assumptions. 


With  this  done,  it  is  possible  to  carry  out  calculations  of  the 
thrust  characteristics  inside  of  ejector  nozzles. 

1.  Concerning  the  common  boundary  between  the  main  flow  and 
the  harmonic  flow(s),  that  is,  the  body  of  revolution  which  is 
formed  by  the  elliptical  curve  is  the  form  of  the  column  of  the 
main  flow. 

2.  During  the  time  in  which  the  flow  inside  the  ejector 
nozzle  is  expanding  and  accelerating,  we  accept  the  fact  that  the 
nozzle  is  in  a  flow  configuration  of  complete  stability,  that  is  to 
say  that  the  working  fluid  does  not  give  rise  to  a  separation  into 
components  and  the  specific  heat  remains  constant. 

3.  There  is  no  heat  exchange  between  the  main  flow  and  har¬ 
monic  flows . 

4.  During  flow,  the  most  important  consideration  is  shock 
wave  losses. 

5.  We  recognize  the  fact  that  the  gas  flow  in  the  ejector  nozzle 
functions  as  an  adiabatic  energy-insulated  flow. 

6.  The  harmonic  flow(s)  act  as  one-dimensional,  isoentropic  flows. 

Figure  1  is  the  shadow  photograph  of  the  flow  spectrum  of  the 
exhaust  of  the  ejector  nozzle  which  we  were  considering  in  three 
different  cases  (these  three  cases  can  be  seen  in  schematic  diagram 
Figure  2).  From  this  set  of  photographs,  it  can  be  seen  that,  under' 
differing  types  of  pressure  loss  ratios,  the  main  flow  always  main¬ 
tains  within  itself  the  existence  of  a  shock  wave  and  following  a 
reduction  in  the  ratio  of  pressure  drop,  the  shock  wave  continously 
moves  in  toward  the  nozzle  sleeve  of  the  ejector  nozzle.  In  this 
condition,  the  shock  wave  also  gains  in  strength.  In  order  to  make 
the  analysis  more  convenient,  we  took  the  flow  configurations  inside 
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(*)  ”  2.56 

(b)  irt-4.7 

(c)  U.7 

oi  —  5% 
oi  —  5% 

£u  —  5% 

Figure  1.  The  flow  configuration  when  different  pressure  drop 
ratios  are  in  effect  for  expansion  type  outer  tube 

sleeves 

ejector  nozzles  for  different  pressure  drop  ratios  and  simplified 
them  by  using  the  models  which  are  shown  in  Figure  2. 

When  the  pressure  drop  ratio  for  the  main  flow  is  larger  than 
the  design  value,  that  is,  when  (*L>nUi)  ,  then,  the  ejector  nozzle 
is  put  into  an  incomplete  expansion  operational  configuration,  and 
the  main  flow  acts  as  an  isoentropic  flow  inside  the  tube  sleeve. 

The  3hock  wave  in  the  main  flow  makes  its  appearance  in  the  free 
Jet  in  the  ejector  nozzle.  Concerning  the  position  of  the  largest 
horizontal  tangential  surface  in  the  column  of  the  main  flow,  this 
may  easily  appear  at  a  certain  place  outside  the  exhaust  of  the 
ejector  nozzle  and  the  actual  position  of  it  is  related  to  the  press¬ 
ure  drop  ratio. 

In  order  to  improve  the  performance  of  ejector  nozzles  under 
conditions  in  which  the  pressure  drop  ratios  are  low,  it  is  only 
necessary  to  have,  when  one  approaches  the  design  pressure  drop 
ratio,  the  main  flow  in  the  ejector  nozzle,  and  it  will  approach 
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complete  expansion.  In  this  type  of  flow  configuration,  the  larg¬ 
est  horizontal  tangential  surface  in  the  main  flow  column  is 
located  on  a  tangential  surface  close  to  the  exhaust  of  the  tube 
sleeve . 

It  makes  no  difference  if  the  main  flow  is  placed  into  an 
incomplete  expansion  configuration  or  into  a  complete  expansion  con¬ 
figuration,  the  tangential  surface  to  the  exhaust  of  the  ejector 
nozzle  is  filled  up  with  lower  harmonic  flow  and  with  supersonic  main 
flow  which  has  not  passed  through  the  shock  wave. 

When  the  main  flow  is  in  the  condition  <  *u,  ,  and  the  main 

jet  nozzle  is  still  in  operation  with  the  pressure  drop  at  levels 
beyond  the  critical  stage,  that  is  ,  then  in  such  a  situa¬ 

tion,  the  ejector  nozzle  is  placed  into  an  over-expansion  configura¬ 
tion.  In  this  type  of  condition,  the  main  flow  leaves  the  main  jet 
nozzle  and  immediately  wraps  itself  around  the  edges  of  the  exhaust, 
that  is,  it  immediately  goes  into  Prandtl-Meyer  expansion.  From  the 
critical  pressure  value  of  the  exhaust  of  the  main  jet  nozzle,  it 
immediately  expands  to  the  point  where  the  edges  of  the  main  Jet 
nozzle  are  placed  in  a  static  pressure  of  harmonic  flow.  After  the 
main  flow  crosses  the  Class  II  expansion  characteristic  curve,  and 
then,  also  crosses  the  Class  I  expansion  characteristic  curve  which 
appears  from  the  opposite  side,  it  is  then,  and  only  then,  that  it 
stops  progressively  expanding.  Because  of  the  fact  that  the  form  of 
the  common  boundary  between  the  main  flow  and  the  harmonic  flow(sl 
is  most  certainly  not  congruent  with  the  surface  form  which  is  re¬ 
quired  by  the  ideal  expansion  of  the  main  flow.  In  this  way,  when 
the  expansion  characteristic  curve  reaches  the  common  boundary  sur¬ 
face,  it  will  reflect  from  the  common  boundary  surface  weak  shock 
waves.  These  weak  shock  waves  gather  together  and  they  form,  at  a 
place  which  is  not  far  from  the  common  boundary  surface,  a  relatively 
strong  shock  wave.  When  the  gas  flow  passes  across  the  shock  wave 
we  just  mentioned,  there  is  created  a  turning  in  the  opposite  direct¬ 
ion.  Moreover,  on  the  largest  tangential  surface  to  the  column  of 
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the  main  flow,  the  gas  flow  also  turns  back  to  the  same  direction 
that  the  axis  has.  It  can  be  seen  from  this  that,  in  configura¬ 
tions  in  which  we  have  excessive  expansion,  the  position  of  the 
largest  horizontal  tangential  surface  to  the  column  of  the  main 
flow  is  placed  inside  the  tube  sleeve  of  the  ejector  nozzle.  More¬ 
over,  on  the  tangential  surface  to  the  exhaust  of  the  ejector  nozzle, 
we  take  the  main  flow  and  divide  it  into  two  parts,  that  is  to  say, 
the  center  flow  which  has  not  passed  through  the  shock  wave  and  the 
exterior  flow  which  has  passed  through  the  shock  wave. 

On  the  basis  of  what  has  been  said  above,  it  is  possible  to  set 
up  the  necessary  equations.  We  recognize  the  fact  that  the  common 
boundary  surface  between  the  main  flow  and  the  harmonic  flow  (.si,  that 
is  to  say,  the  main  flow  jet  column  is  related  to  the  body  of  revo- 
lution  which  is  formed  by  the  elliptical  curve  and  its  gear  profile 
as  well  as  to  the  pressure  drop  ratio  and  the  flow  ratio  ui. 

Because  of  these  facts,  the  cylinder  of  the  main  flow  jet  can 
be  described  by  the  use  of  this  equation 

- — _ +.  JtL  -  i  Cl 1 

(R\.4  -  R\)Ri.J  'g1  dR\  R\ti 

In  this  equation,  Rzzd  is  the  radius  of  the  largest  horizontal 
tangential  surface  to  the  column  of  the  main  flow  and,  on  the  basis 
of  experimental  data,  it  is  possible  to  adjust  it  to  be 

_  rp ) 

R„4  ““  RmV  1  +  K(*L  —  *LL/) 


Experimental  results  demonstrate  that,  within  the  range 
«L <  g  —  io,  uj  “  o.oi -- u.04  it  is  possible  to  select  K  =  0.18. 

The  angle  0  in  Equation  (.1)  is  the  angle  of  slope  of  the  cylin¬ 
der  of  the  main  jet  at  the  exhaust  of  the  main  jet  nozzle.  This  is 
determined  by  the  expansion  of  the  pressure  which  is  found  on  the 
tangential  surface  of  the  main  flow  at  the  exhaust  of  the  main  jet 
nozzle  expanding  into  the  pressure  which  is  found  in  the  surrounding 
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Figure  2.  Simplified  models  of  ejector  nozzles  under  different 
flow  configurations. 

1 — harmonic  flow;  2 — main  flow;  3 — excessive  extension;  4 — harmonic 
flow;  5 — cylindrical  surface  of  the  main  jet;  6 — shock  wave;  7 — 
complete  expansion;  8 — incomplete  expansion 


harmonic  flows,  that  is  to  say,  d  —  f(p,/p ?)  (when  p,<+p?  ,  it  is 

it  is  possible  to  select  9  —  fcpf/p!)).  Due  to  the  fact  that  the 
main  flow,  at  the  exhaust  of  the  main  jet  nozzle,  suddenly  expands 
around  the  outside  obtuse  angles,  the  main  flow,  which  expands  in 
three  dimensions,  can  be  usefully  illustrated  in  two  dimensions. 


From  this  it  is  possible  to  fix  the  distance  x1  which  separates 
the  largest  horizontal  tangential  surface  of  the  main  jet  cylinder 
from  the  edge  of  the  exhaust  of  the  main  jet  nozzle. 

Ki)/'g*R.  C3) 
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In  an  incomplete  expansion  type  situation,  x-^  >  L,  that  is  to 
say  that  the  largest  tangential  cross  section  of  the  cylinder  of  the 
main  flow  is  located  outside  of  the  exhaust  of  the  ejector  nozzle. 
The  radius  of  the  main  flow  cylinder  on  the  tangential  plane  of 
the  exhaust  of  the  ejector  nozzle  is 


(xt-  LytK>6R\ 

-  R\)R\.d 


(.4) 


Once  again,  on  the  basis  of  the  continuity  conditions  of  the 
main  flow  contained  in  the  relationship  tt(X),  it  is  possible  to 
obtain  all  the  physical  parameters  as  well  as  the  tangent  planes  of 
the  path  of  flow  on  the  tangent  plane  of  the  exhaust  of  the  ejector 
nozzle  for  the  main  flow. 


In  this  sort  of  operational  configuration,  the  tangent  plane 
which  is  occupied  by  the  harmonic  flow  on  the  tangent  plane  to  the 
exhaust  of  the  ejector  nozzle  is 

A„  -  *(/«  -  R‘„)  J 


Concerning  the  flow  speed  of  the  harmonic  flow  at  the  exhaust 
of  the  ejector  nozzle,  this  can  be  obtained  from  the  equation  for 
the  amount  of  flow 


*<*»> 


P?  Altm,C , 


(6) 


In  a  tube  sleeve  with  a  contracted  form,  the  greatest  speed  of 
flow  of  the  harmonic  flow  at  the  exhaust  is  —  l.  Generally 

speaking,  when  the  pressure  drop  ratio  and  the  ratio  of  the  amounts 
of  flow  are  not  too  large,  that  is 

(«.  -  ( G,/G.)y/TTIT •) 

then,  the  flow  speed  of  the  harmonic  flow  in  the  tangent  plane  to 
the  exhaust  of  the  tube  sleeve  is  always  subsonic. 
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In  this  type  of  case,  that  is,  when  the  main  flow  is  in  a 
condition  of  complete  or  incomplete  expansion,  and  the  flow  speed 
of  the  harmonic  flow  at  the  exhaust  has  not  reached  a  critical 
value,  the  actual  thrust  which  Is  produced  at  the  ejector  nozzle 
is  P  =  F_  +  F  ,  that  is  to  say 

F  —  —  n*»)  —  pn  1  +  —  pul  ^ 


In  this  equation,  n  is  the  corrected  thrust  parameter  as  de¬ 
rived  from  the  ideal  theoretical  parameter  which  does  not  take  into 
consideration  the  effects  of  boundary  layers.  For  the  respective 
cases  of  the  main  and  harmonic  flow,  this  parameter  is 


y»  ** 


c a  *  c 4g 

c,  •  c* 


(8) 


Cd,  C  are  respectively  the  parameters  for  the  amount  of  flow 
and  momentum  losses  as  determined  by  experimentation.  When  the 
main  flow  is  in  an  excessive  expansion  configuration,  as  a  function 
of  a  drop  in  the  pressure  drop  ratio,  the  shock  wave(s)  moves  up  the 
path  of  the  gas  flow  as  far  as  the  upper  reaches  of  the  flow. 

Because  of  this,  in  a  situation  in  which  there  is  a  relatively  low 
pressure  drop  ratio,  the  shock  wave  Is  capable  of  entering  as  far 
as  the  inside  of  the  tube  sleeve  of  the  ejector  nozzle.  At  such  a 
time,  the  largest  horizontal  tangent  plane  to  the  cylinder  of  the 
main  flow  i3  located  within  the  tube  sleeve.  The  shock  wave  takes 
the  main  flow  and  divides  It  up  into  two  parts  on  the  tangent  plane 
to  the  exhaust  of  the  tube  sleeve,  that  is  to  say,  into  the  center 
flow  in  front  of  the  shock  wave  and  the  external  flow  behind  it. 


The  flow  speed  of  the  central  flow  in  front  of  the  shock  wave 
is  obtained  from  the  continuity  condition 

*(*„)  -(*./*., JJ  (51 


Its  static  pressure  is 


U01 
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X '  for  the  main  flow  behind  the  shock  wave  is  determined  by 

Z  0 

the  strength  of  the  shock  wave.  The  reason  for  this  is  that,  when 
one  is  dealing  with  the  type  of  situation,  one  can  recognize  the 
fact  that  after  this  section  of  the  flow  goes  through  the  cone  of 
the  shock  wave,  its  static  pressure  goes  up  to  that  of  the  atmos¬ 
pheric  pressure  PH  on  the  external  boundary.  The  shock  wave  angle 
0g  can  be  determined  from  the  pressure  rise  ratio  ph/pzc  and  the 
speed  in  front  of  the  wave  X  .  Because  of  this,  it  is  possible  to 
make  use  of  the  related  equations  which  are  already  known  or  the 

f 

graphs  of  these  in  order  to  figure  out  0s,  Xzc  as  well  as  the  para¬ 
meter  of  total  pressure  losses  a. 


In  this  type  of  operational  configuration,  the  radius  of  the 
main  flow  cylinder  on  the  tangent  plane  at  the  exhaust  of  the 
ejector  nozzle  is 


K.' 
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By  making  use  of  geometrical  relationships,  we  can  figure  that 
the  area  of  the  tangent  plane  which  Is  occupied  by  that  part  of 
the  main  flow  which  has  not  gone  through  the  shock  wave  is 

Rzt\~Rt*i — (i.- C.H) 

■  i»«i  “  *R\t i 


The  area  of  the  tangent  plane  at  the  exhaust  of  the  tube  sleeve 
which  Is  occupied  by  that  part  of  the  main  flow  which  has  passed 
through  to  the  other  side  of  the  shock  wave  Is 

-V.„)  Cl2i 


Because  of  these  facts,  when  one  Is  operating  in  a  situation  of 
excessive  expansion,  the  actual  thrust  which  is  produced  by  the 
ejector  nozzle  is 


F  -  P ...  +  F,m  +  F, 


that  is  to  say. 


F  “  i)  +  +  ^7~T  Vu  ri 

'•(*,.)  113  2 

—  Ph(‘4  >f i  ^  'til  i  ^7 if  I  ^  w^i# } 

In  this  way,  when  one  considers  ejector  nozzles  with  already 
given  geometrical  forms  to  their  external  shapes  as  well  as 
/>?•.  p?„  77,  77,  u>  the  main  and  harmonic  flows  and  the  properties  of 
the  main  and  harmonic  flows  k,  m,  then  it  is  possible  to  calculate 
the  actual  thrust  which  is  produced  when  the  ejector  nozzle  under 
consideration  is  operating  under  different  configurations.  And  from 
this,  we  can  go  on  to  figure  out  the  relationship  between  this 
thrust  and  the  pressure  drop  ratio,  that  is  to  say,  the  thrust  per¬ 
formance  inside  the  ejector  nozzle  under  consideration. 

(Two)  The  actual  thrust  inside  of  the  ejector  nozzle 

If  one  looks  at  Figure  3,  it  can  be  seen  that  the  solid  wall 
surfaces  which  the  main  and  harmonic  flow  wrap  themselves  around  are 
all  show  the  presence  of  the  formation  of  boundary  layers,  Concern¬ 
ing  the  area  of  mixing,  from  the  trailing  edges  of  the  main  jet 
tube  onward,  when  the  main  flow  and  the  harmonic  flow  begin  to  come 
in  contact  with  each  other,  they  flow  on  smoothly  in  the  direction 
of  the  air  flow,  expanding  as  they  go  until  they  reach  the  exhaust 

of  the  ejector  nozzle.  The  thickness  of  the  area  of  mixing  is  6 

**  c 
and  its  thickness  of  momentum  loss  Is  SQ  .  If  we  make  use  of  the 

laws  governing  momentum,  we  accept  the  condition  of  ignoring 
the  shear  forces  on  the  boundaries  of  the  area  of  mixing.  Concern¬ 
ing  the  flow  within  the  area  of  mixing.  It  is  possible  to  use  this 
type  of  impulse  equation  in  order  to  describe  the  situation 

cm 

dd”/B“dx  +  dp./p.dx  +  (2  +  «•/«**)  •  dVjV,dx  -  0. 

If  we  take  the  parameter  for  the  form  of  the  surface  H  -  <5*/6** 

for  purposes  of  expression  and  we  take  the  fixed  values  of  Hpj  « 

(H  +  H_)/2  for  purposes  of  substitution,  then  the  equation  above 
z  z  c 

can  be  directly  integrated  and,  after  it  has  undergone  simplification 
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we  can  obtain  the  thickness  of  momentum  loss,  , 

V 

plane  at  the  exhaust  of  the  ejector  nozzle  is 


on  the  tangent 


-  p./p? 


(p«/p.*)  • 

-  [2/U,  +  i)3|/*,_l  •  -  -  .  -~~x - 

(  1  - 

V  +  l  ' 


Prom  the  equation  above,  it  can  be  seen  that,  as  far  as  the 
thickness  of  momentum  loss  at  the  exhaust  of  the  ejector  nozzle  goes, 
besides  the  fact  that  it  is  related  to  the  way  in  which  the  main 
flow  expands,  this  thickness  is  also  determined  by  the  parameter  for 
the  average  surface  form  Hpc-  In  order  to  make  a  precise  determina¬ 
tion  of  the  thickness  of  momentum  loss  at  the  exhaust  of  the  Jet 

i* 

nozzle  for  the  region  of  mixing,  that  is  <$  ,  we  assume  that  the 

distributions  of  speed  and  temperature  in  the  area  of  mixing  are 


If  one  makes  use  of  the  relationships  involving  U,  T  in  the 

equations  above  as  well  as  of  p/p_  *  T  /T,  then  it  is  possible  to 

z  z 

obtain 


Figure  3.  The  actual  thrust  in  the  ejector  nozzle 
1 — boundary  layer;  2 — area  of  mixing 
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After  adjustments  are  made,  we  get 

A  1 

1  A 

2  B  V'l  -~Z>2 

In  these  equations 


8 


^-i  +  W-i  +  a,/a.  •  Vrr/r? 
a  - 1  +  r,/r.  -  i  +  tt  •  r(i,)/7.*r(;.,) 

o  -  1  -  7-./y.  _  t  -  r,*ra)/77rU.) 

1  +  T,/T,  1  +  7,*r(A,)/ rfr(  A,) 


Because  of  this,  it  is  possible  to  get  a  precise  value  for 

£  ft 

the  thickness  of  momentum  loss  6C  at  the  exhaust  of  the  Jet  nozzle 
for  the  area  of  mixing,  and  from  this,  one  can  go  on  to  get  a  pre¬ 
cise  value  for  the  parameter  of  adjusted  thrust  n  for  a  situation 
in  which  one  is  taking  mixing  into  consideration. 

3.  A  COMPREHENSIVE  STATEMENT  CONCERNING  EXPERIMENTATION 

COne )  A  short  Introduction  to  experimental  equipment 

One  can  see  a  picture  of  the  test  bed  in  photograph  4  and  a 
schematic  diagram  of  it  in  Figure  5.  It  is  composed  of  such  compo¬ 
nents  as  curved  air-guide  tubes  and  main  flow  passages.  Concerning 
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the  main  Jet  tubes  and  external  tube  sleeves  which  were  used  in 
the  experiments,  both  types  of  components  were  screwed  into  con¬ 
nections  on  the  test  bed.  The  main  flow  passages  and  the  curved 
air-guide  tubes  used  floating  tube  type  connectors.  The  interval 
between  components  was  0.6  mm,  and  this  interval  was  maintained  in 
order  to  make  it  easier  for  the  components  involved  to  float  free 
under  the  effects  of  applied  thrust.  The  harmonic  flow  passages 
and  the  main  flow  flow  paths  were  connected  together,  and  this  was 
done  in  order  to  facilitate  the  process  of  making  the  pressure  in 
the  harmonic  flow  passages  the  same  as  that  in  the  main  flow. 

Concerning  the  dry  fresh  air  which  is  supplied  by  the  air  move¬ 
ment  ducts,  the  largest  part  of  it  is  put  out  through  the  main  Jet 
nozzle.  A  small  part  of  it  enters  the  outside  tube  sleeve  and  forms 
the  harmonic  flow  used  in  the  experiments.  The  amount  of  the  har¬ 
monic  flow  is  changed  by  the  exchange  limiter  nozzle. 

The  test  bed  was  equipped  with  weighted  lever  mechanisms  and 
these  were  used  to  make  standard  measurements'  of  the  system  of  forces 
involved  in  both  static  and  dynamic  states. 

(.Two)  Measurement  parameters  and  instruments 

We  made  use  of  total  pressure  sensors  in  order  to  measure  the 
total  pressures  of  the  main  flow  and  the  dynamic  flow  being  consi¬ 
dered.  Besides  this,  there  were  drilled,  in  the  wall  of  the  harmonic 
flow  ring  cavity  and  in  the  wall  surface  of  the  air-guide  tubes 
below  the  limiter  nozzle,  static  pressure  holes  1  mm  in  diameter, 
and  these  holes  were  used  for  making  precise  measurements  of  the 
static  pressure  of  the  harmonic  flow. 

All  the  measurements  of  pressure  went  through  a  pressure  sensor 
(type  SYD),  circulating  measurement  devices  CXJ-1QQ1  and  the  digital 
recorder  CLS-5).  It  was  then  directly  displayed  as  numerical  press¬ 
ure  values  or  printed  out  as  pressure  data. 


Use  was  also  made  of  force  measuring  sensors  (BLR-1),  carrier 
wave  amplifiers  (LR-1)  and  integral  type  numerical  voltage  tables 
(DS-13)  as  well  as  of  printers  (SY-3)  in  order  to  make  measurements 
of  the  thrust  affecting  the  test  components. 

(Three)  Test  components  and  experimental  sequence 

We  made  tests  of  several  score  of  main  jet  tubes  with  differ¬ 
ent  exhaust  diameters  (60-90  mm)  and  different  angles  of  narrowing 
(a  from  10°-14°) 


Figure  4.  The  test  bed 

We  also  made  tests  of  several  score  of  differnt  forms  of  exter¬ 
nal  tube  sleeves.  These  were  various  sleeves  with  different  forms 
of  narrowing,  narrowing  and  then  expanding,  sleeves  with  different 
diameter  ratios  (D  /D  ) ,  different  spacing  ratios  CL/D),  differ- 
ent  angles  of  sleeve  contraction  a  and  sleeve  expansion  8  as  well 
as  other  similar  variable  parameters.  All  these  various  forms  were 
tested. 

All  the  testing  was  done  within  a  range  in  which  the  ratio  of 
the  amounts  of  flow  w  is  between  ll-7%  and  the  pressure  drop  ratio 
is  between  1.5  and  14.  Moreover,  for  each  value  of  the  ratio  of  the 
amounts  of  flow,  continuous  experimentaion  was  carried  out  on  the 
effects  of  changes  in  the  pressure  drop  ratio. 
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4.  A  COMPARISON  OP  THE  CALCULATED  VALUE  OP  THRUST  AND  THE  ' 
EXPERIMENTALLY  MEASURED  VALUE 

We  selected  the  experimental  data  for  tube  sleeves  with  seven 
different  types  of  narrowing  forms  and  with  a  spacing  ratio  L/D 

z 

from  0.58-0.6,  and  we  made  a  comparison  between  these  results  and 
those  thrust  values  which  were  obtained  by  calculation  using  the 
methods  put  forward  in  this  article.  The  results  of  this  compari¬ 
son  are  shown  in  Figure  6.  It  can  be  seen  from  this  illustration 
that  within  a  very  large  range  of  values  for  the  pressure  drop  ratio 
of  the  main  flow,  these  two  sets  of  values  agree  with  each  other 
very  well.  Figure  7  shows  the  error  between  the  calculated  values 
and  the  experimentally  determined  values  as  well  as  the  relation¬ 
ship  to  the  pressure  drop  ratio.  It  can  be  seen  from  this  illustra¬ 
tion  that  in  practical  applications  the  greatest  encountered  any 
where  is  not  larger  than  4£  and  the  distribution  of  this  error  is 
relatively  concentrated.  It  is  only  in  the  case  in  which  the  ejec¬ 
tor  nozzle  is  placed  in  an  excessive  operational  configuration  that, 
due  to  the  fact  that  the  form  of  the  flow  inside  of  the  tube  sleeve 
is  very  complicated,  there  is  caused  a  significant  reduction  in  the 
accuracy  of  estimates.  Because  of  the  fact  that  this  article  consi¬ 
ders  the  characteristics  of  the  actual  flows  inside  of  the  ejector 
nozzles  under  different  operational  configurations,  it  raises  the 
accuracy  of  the  estimates  of  thrust  over  the  entire  range  of  pressure 
drop  ratios.  Because  of  this  fact,  this  method  is  capable  of  func¬ 
tioning  as  something  to  which  engineers  can  refer  when  selecting  the 
initial  design  plans  for  the  structure  of  ejector  nozzles, 

5 .  CONCLUSIONS 

The  method  of  analysis  of  the  characteristics  of  the  inside  of 
ejector  nozzles  which  is  put  forward  by  this  article  as  well  as  the 
corresponding  methods  for  estimating  the  values  of  these  character¬ 
istics  have  been  verified  by  going  through  preliminary  testing. 

These  methods  are,  within  a  very  wide  range  of  pressure  drop  ratios 
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II  *3 


Figure  5.  A  simplified  diagram  of  the  test  bed 

1 — gas  leakage;  2 — force  measuring  sensor;  3 — flow  limiter  jet 
nozzle;  4 — harmonic  flow  collection  chamber;  5 — wall  surface 
pressure  measurement;  6 — main  Jet  nozzle;  7 — flow  rectifier  grid; 

8 — harmonic  flow;  9 — pressure  measuring  rake;  10 — intake;  11 — 
curved  intake  tube;  12 — free  flow;  13 — harmonic  flow  static  press¬ 
ure;  14 — external  tube  sleeve;  15 — strength  measurement  calibration 
structure;  16 — forward  sprung  plate;  17 — backward  sprung  plate 


ttl  (approximately  from  2.5-14)  and  ratios  of  amounts  of  flow 
(approximately  capable  of  estimating  calculated  values 

of  thrust  which  are  in  satisfactory  agreement  with  the  experiment¬ 
ally  measured  values  for  the  same  quantities. 


By  making  use  of  these  same  types  of  methods,  it  is  also  poss¬ 
ible  to  deduce,  for  different  operational  configurations  of  ejector 
nozzles,  the  situation  concerning  changes  in  the  structure  of  the 
shock  wave  in  the  main  flow.  This  helps  us  to  understand  the  nature 
of  the  flow  in  the  ejector  nozzles  being  studied. 
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Figure  6.  A  comparison  between  the  calculated  values  for 
ejector  nozzle  thrust  and  the  experimentally  measured  values 
corresponding  to  them. 


Figure  7.  The  errors  between  the  calculated  values  for  thrust  in 
the  ejector  nozzles  being  studied  and  the  experimentally  measured 
values  for  the  same  quantity. 


Figure  7  notes: 


1/0,  -U.5W  a>  —  2.5¥%  O 
U  —  3.5%  A 
<ii  —  5 .0%  ® 

L/D,  -  0.6033  m  -  1.5%  A 
u-2.5%  X 
oi-3.5%  □ 
LID,  -0.4676  ui-  3.5%  » 
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ANALYTICAL  AND  EXPERIMENTAL  INVESTIGATION  OF 
SUPERSONIC  EJECTOR  NOZZLE 


Wu  Da  Yang  Zai-ming  Zhang  Rung 
(Beijing  Institute  of  Aeronautics  and  Astronautio*) 


Abstract  *r 

Tit  this  paper  we  introduce  some  results  of  test  of  ejector  nozzle,  and  propose  in 
compliance  with  these  results  a  model  of  the  flow  in  ejector  nozzle,  and  so  we  establish 
ait  approximate  method  to  count  the  thrust  characteristics  of  it.  At  last  we  compare  the 
results  with  data  of  test  and  discover  that  they  are  good  in  agreement. 


36 


B1  TYPE  TRANSONIC  AXIAL-PLOW-COMPRESSOR  FIRST  STAGE  ROTOR 
BLADE  DESIGN,  TEST  AND  APPLICATION* 

Zhou  Wei-guo 

(.Shenyang  Aeroengine  Company) 

SUMMARY 

This  article  describes  the  process  of  experimentati  .1  and 
redesign  which  was  carried  out  in  order  to  eliminate  failures  due 
to  cracking  of  I  stage  type  engine  blades.  This  article  is  com¬ 
posed  altogether  of  three  sections.  The  first  section  is  one  which 
deals  with  the  methods  used  to  redesign  the  blades  in  question. 

The  second  section  contains  the  results  of  experimentation.  And 
the  third  section  is  an  appraisal  of  the  application  of  these  blades. 

FOREWORD 

The  B1,  from  the  time  of  its  testproduction ,  both  inside  the 
factory  and  outside  the  factory,  developed  many  instances  of  fail¬ 
ures  due  to  cracking  and  breaking  of  blades  in  the  first  stage  oper¬ 
ations  of  the  engine.  In  the  factory,  there  was  one  case  of  bending, 
four  cases  of  cracking.  Out  in  the  field,  there  was  one  case  of  a 
major  accident  and  five  cases  of  cracking  of  these  components,  which 
severely  influenced  the  safety  of  aircraft  in  the  Peoples  Air  Force. 

At  that  time,  our  analysis  was  this:  We  thought  that  the  fifth 
batch  of  engines  had  common  operational  curves  which  were  relatively 
close  to  the  surge  line  and  the  surge  margin  was  small.  Because  of 
this,  we  carried  out  a  series  of  dynamic  and  static  frequency  tests. 
With  the  approval  of  the  factory,  we  produced  a  fifth  batch  of 
thicker  blades  (3imply  called  fifth  batch  thick  blades).  When  dyna¬ 
mic  and  static  frequency  tests  were  made  to  determine  these  fre¬ 
quencies,  the  results  of  the  experiments  showed  no  obvious  differ¬ 
ences  between  the  basic  natures  of  the  normal  fifth  batch  blades  and 
the  thicker  fifth  batch  blades. 


This  article  was  read  at  the  2nd  All  China  Engineering  Thermo- 
physics  Technology  Conference  in  Hangzhou  in  November  of  1978, 
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In  1965,  we  made  experimental  tests  for  a  sixth  batch  of 
engines.  The  first  stage  operational  blades  and  the  flow  recti¬ 
fier  blades  in  the  fourth  and  fifth  stages  saw  some  relatively  large 
changes  from  those  of  the  previous  model.  Among  these  changes  were 
the  facts  that  the  operational  blades  in  the  first  stage  compress¬ 
ors  Ccalled  "first  flow"  for  short)  had  their  arc  lengths  ex¬ 
tended  by  21-26  mm.  The  number  of  blades  was  changed  from  31  to 
24.  The  arc  length  of  the  fourth  stage  rectifier  blades  was  re¬ 
duced  by  3-4  mm.  And,  the  maximum  thickness  of  the  fifth  stage 
rectifier  blades  was  increased  0.15-0.24  mm.  The  flow  parameters 
resulting  from  these  three  changes,  with  the  exception  of  the  five 
cross  section  of  the  first  stage  compressor  blade  tips  pXKj  plK  all 
exhibited  changes.  All  the  other  parameters  were  the  same  as  they 
were  with  the  old  form.  In  order  to  completely  eliminate  malfunc¬ 
tions  due  to  problems  with  the  first  stage  compressor  components, 
in  September  of  1965,  work  was  begun  on  redesigning  the  first  stage 
compressors.  In  1970  and  in  the  years  1975-1976,  we  made  more  than 
20  tests  and  these  experiments  are  introduced  below. 

1.  DESIGN  OF  THE  OPERATIONAL  BLADES  OF  THE  FIRST  STAGE  COMPRESSORS 
(CONFIGURATION) 

1.  The  arc  length  of  the  first  stage  compressor  blades  in  the 
sixth  batch  of  engines  was  lengthened  by  21-26  mm  over  the  length 
as  it  was  in  the  fifth  batch  of  engines.  The  number  of  blades  was 
reduced  from  31  to  24.  The  blade  configuration  parameter  b/t  was 
the  same  as  it  was  with  the  fifth  batch.  And  the  cross  sections 
Bik>  82k  of  the  five  blade  tips  were  larger  than  those  for  the  fifth 
batch  of  engines  by  lo56,-lo07,.  The  results  of  these  experiments 
were  these.  The  temperature  of  the  exhaust  gases  of  the  sixth  batch 
of  engines  as  well  as  their  rates  of  fuel  consumption  were  relatively 
higher  than  the  corresponding  values  for  the  fifth  batch  of  engines. 
From  this  fact,  it  can  be  seen  that  the  aerodynamic  performance  of 
the  fifth  batch  of  engines  was  superior  to  that  of  the  sixth  batch 
of  engines.  When  we  gave  consideration  to  the  cracking  malfunctions 
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associated  with  the  first  stage  compressor  blades  of  the  fifth  batch 
of  engines,  we  referred  to  the  arc  length  and  maximum  thickness  of 
the  blades  of  the  first  stage  compressors  of  the  sixth  batch  of 
engines  in  order  to  select  a  first  stage  compressor  blade  config¬ 
uration  designed  on  the  basis  of  these  data  and  the  form  which  was 
designed  was  the  1-02  blade. 

This  designing  process  was  carried  out  in  the  years  1965-66. 

The  calculations  of  intake  flow  fields  which  are  contained  in  this 
report  were  done  on  the  basis  of  the  two-dimensional  type  calcula¬ 
tion  method.  The  blade  parameters  used  can  be  seen  set  out  in  the 
table  below. 

One  can  see  from  the  data  in  the  table  below  that  the  fact  that 

the  arc  length  is  increased  and  the  Reynolds  number  is  increased 

broadens  the  range  of  operational  stability  of  the  engine.  The 

value  (T  „  for  the  sixth  batch  of  engines  is  smaller  than  that  for 
max 

the  fifth  batch,  this  is  an  advantage  for  the  supersonic  and  tran¬ 
sonic  blade  performance  involved.  A  cross  section  of  the  I  blade  of 
the  sixth  batch  has  a  level  straight  section  to  it.  Below  the  II 
cross  section  the  intake  has  a  shorter  straight  section  to  it.  The 
blade  trough  has  a  section  which  gradually  widens  and  has  no  throat 
section  to  it.  This  sort  of  configuration  for  supersonic  blades 
does  not  correspond  to  what  is  normally  dope  with  this  type  of 
design  and,  in  the  design  process,  adjustments  are  needed. 

2.  Operational  configurations  of  first  stage  blades  (.see 
Table  2). 

3.  The  method  of  configuring  the  shape  of  supersonic  blades 
(appropriate  for  use  with  values  of  M  *  1.15-1. ^01. 

(1)  Draw  a  straight  line  OQ  which,  with  the  quota  lines,  forms 
01K  ^-see  PiSure  D*  In  that  figure  OQ  —  *(l— iMr'/A)  and,  in  this 
relationship 
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TABLE  1.  A  table  of  the  parameters  for  the  blade  configurations 
for  the  fifth  and  sixth  batches  of  engines 


s  section 

II 

I 

11 

1U 

IV 

V 

VI 

VII 

V1U 

IX 

X 

XI 

H 

110.5 

105.1 

279.7 

254.1 

228.8 

201.5 

178.05 

165.4 

n .  6*/ 

80.53 

85.66 

01.47 

81.05 

78.76 

76.  r 

73.61 

73.02 

b 

116 

114 

110 

107 

101.8 

101.4 

98.1 

94.8 

91.5 

(92.6) 

(92) 

(91) 

Hi 

1. 107 

1.322 

1.383 

1.472 

1.572 

1.70 

1.851 

2.04 

2.18 

29*20' 

31*23' 

16*57' 

43*45' 

51*49' 

62*9' 

76*18' 

96*6' 

120“ 

0*IC 

12*10' 

18*30' 

45* 

53*45' 

62“ 

78* 

97*30' 

115* 

(131*30') 

(150*45') 

— 

25*15' 

27*14' 

31*26' 

15*6' 

17*28' 

42*55' 

47*7' 

51*51' 

54*10' 

57*22' 

61*42' 

65*16' 

0i* 

29* 

33* 

36*10' 

39* 

44* 

47*15' 

51*10' 

54*15' 

58*45' 

63* 

2.75 

2.93 

3.72 

4.46 

5.23 

6.07 

6.97 

7.90 

8.23 

'■* 

2.19 

3.03 

3.80 

4.63 

5.30 

6.19 

7.15 

7.70 

(8.37) 

(8.70) 

•*c% 

64.6 

65.2 

56.2 

53.8 

53.1 

54.3 

56.8 

57.6 

61 

55.7 

54.7 

51.5 

57.7 

55.3 

57.1 

57 

62 

65 

- 

5*45' 

6*12' 

7*15' 

7*27' 

7*16' 

7*06' 

7*30' 

8*18' 

•« 

3*30' 

3* 

4* 

4* 

5*15' 

5* 

6*15' 

4*33' 

4*33' 

5* 

5*42' 

7*42' 

9*54' 

,2*6' 

13*54' 

•» 

3*30' 

4*45' 

5*30' 

7*30' 

9* 

11* 

11*15' 

0.34 

0.35 

0.37 

0.40 

0.43 

0.46 

0.50 

0.57 

0.63 

rl 

0.35 

0.37 

0.40 

0.43 

0.46 

0.50 

0.57 

0.63 

0.55 

0.56 

0.58 

0.62 

0.66 

0.69 

0.72 

0.79 

0.90 

rl 

0.40 

0.45 

0.50 

0.54 

0.55 

0.60 

0.85 

0.90 

NOTE:  The  parameter  values  given  on  top  are  those  for  the  fifth 
batch  and  the  values  given  underneath  are  for  the  sixth 
batch  of  engines. 

(2)  If  one  draws  a  straight  line  Q02  past 
forms,  with  the  rating  line  involved,  the  angle 
the  next  step. 

(3)  Draw  the  tangent  circle  to  the  leading  edge  quota  line  and 
the  line  OQ.  The  radius  of  this  circle  is  r1 . 

(4)  Go  through  the  center  of  the  forward  tangent  circle  0^, 
Take  b — (.r^  +  r2)  to  be  the  radius  at  Q02  and  get  the  point  02< 


the  point  Q  which 
S2K,  then  go  on  to 
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TABLE  2.  Original  data  for  the  1-02  blade  form 


0 

1 

n 

m 

IV 

V 

VI 

vn 

VQI 

R 

339 

330.5 

305.1 

279.7 

254.3 

228.8 

203.3 

178.03 

165.4 

t 

116 

114.38 

110.61 

107.79 

104.66 

101.87 

98.61 

95.12 

94.4 

t\K 

25*13' 

27*14' 

31*26' 

33*6' 

37*28' 

42*55' 

47*7' 

51*51' 

54*30' 

fi.K 

29*20' 

31*23' 

36*57' 

43*43' 

51*49' 

62*09' 

76*1J' 

96*6' 

120* 

Cmm 

2.59 

2.59 

3.45 

4.15 

4.95 

5.50 

6.25 

6.90 

7.20 

Xc% 

64.6 

65.2 

56.2 

53.8 

53.1 

54.3 

56.8 

57.6 

57 

•i 

2*24' 

2*42' 

3*39' 

4*30' 

5*06' 

5*30' 

5*45' 

6* 

6*6' 

If 

5*42' 

3*48' 

4*21' 

5*lg' 

6*54' 

8*42' 

8*54' 

9* 

9* 

It 

53 

49 

'  37 

25.45 

0 

1 

0 

0 

0 

0 

b 

106.3 

UT 

110 

22.75 

0 

0 

0 

0 

0 

hi 

0.49 

0.59 

0.55 

0.62 

0.52 

0.67 

0.61 

0.59 

0.59 

0.53 

0.48 

0.65 

0.59 

0.54 

- 

0.34 

0.35 

0.37 

0.40 

0.43 

0.46 

0.50 

0.57 

0.63 

0.40 

0.40 

0.43 

0.50 

0.54 

0.55 

0.60 

0.85 

0.90 

1.36 

1.33 

1.23 

1.18 

1.09 

1.04 

0.98 

0.91 

0.87 

Key:  1 — back;  2 — dish;  3 — two  straight  lines;  4 — two  straight  lines; 
5 — one  straight  line 


C5)  Going  through  point  02  draw  a  circle  with  radius  r2  ■ 

C6)  Join  0202  the  value  0^  select (xfi — r-^to  get  the 

point  M. 


(.7)  Draw  e^,  e2.  is  formed  by  drawing  a  line  form  OQ  to 

QA2.  e2  is  the  supplementary  line  segment  from  Q02- 

(8)  Extend  OQ  back  away  from  the  blade  and  e2  will  connect 
with  point  B-^  on  one  side  of  angle  e^.  Moreover,  on  OQ  chose 
1.  *  0A* 


C 9)  From  A^  B^,  form  a  triangle  and  select  F  Cconsult  Table 
3)  as  a  secondary  curve. 
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Figure  1.  A  schematic  diagram  of  a  supersonic  blade  configuration 

CIO)  At  M  draw  a  perpendicular  O^Og.  On  this  line,  draw  a 
circle  c nax  tanSent  to  the  curve  of  the  back  of  the  blade. 

Ql)  Extend  in  the  direction  of  the  dish  of  the  blade 

and  select  UA0  *  1  . 

2  P 

(12)  Extend  e2  in  the  direction  of  the  dish  of  the  blade  so 

that  line  segment  PC2  meets  with  UA2  at  point  B2 . 

Draw  a  secondary  curve  in  the  triangle  formed  from 

Check  out  the  location  of  the  point  at  which  the  thick¬ 
ness  is  largest.  If  the  position  of  has  moved  beyond  the 

designed  value  for  it  by  5%,  then  select  new  values  for  ,  e2. 
Begin  the  recalculation  process  from  procedure  (.7)  and  continue  it 
until  one  gets  a  satisfactory  position  for  X  within  an  acceptable 
tolerance  C5X  ?c) 

4.  The  design  of  transonic  blade  forms  (.secondary  curve 
configurations ) . 

(1)  On  the  basis  of  experience,  select  the  position  of  great¬ 
est  curvature  of  the  center  line  Tj.,  At  the  various  cross  sections 
of  the  type  engines  ■  0.5,  F  ■  0.57. 
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(.13) 


A2®2^2  * 


(14) 


p  +  |  —  2X, 

(2)  On  the  basis  of  the  formula  tg<p’  lg  p‘  p  —  i  +  2X ,  as  wel1 
as  0  "  <?i  +  fi  In  combination  it  is  possible  to  get  values  for 
the  angles  ?i>  <>»« 

(.3)  r  —  ^,*4-9,  (according  to  the  angle  of  installation 
angle ) . 

(4)  From  b— r,~ r„  <p„  <p,  draw  the  triangle  ABC  and  check 

the  table  to  find  the  center  line. 

(5)  At  points  A,  C  make  the  circles  r1,  r2  and  the  acute  angles 
£  1 ,  e2<  Then  solve  for  the  tangent  points  A^,  A2  C-^,  C2  between 

the  two  sides  of  the  acute  angle  and  the  circles  r^,  r2 .  After 
this,  extend  the  two  sides  of  the  angle. 

(6)  On  the  centerline,  make  a  perpendicular  to  AC  passing 
through  point  .  This  line  intersects  the  centerline  at  point  M. 

V 

If  one  draws  a  diameter  at  point  M,  then  it  defines  the  circle  . 

(7)  One  then  draws  a  secondary  curve  through  the  triangles 
A1®1C1  and  A2B2C2  whlcil  ls  tangent  to  the  circle  C^^  (_F  is  selected). 

5.  Calculation  of  the  intake  area  flow  field  cross  sections 
in  supersonic  applications. 

Calculations  of  the  intake  area  flow  field  make  use  of  a  posi¬ 
tive  shock  wave  model.  These  calculations  are  carried  out  on  the 
basis  of  the  methods  outlined  In  References  [13  and  [2]  and  will 
not  be  discussed  in  this  article. 

2.  OPERATIONAL  BLADE  TESTS  ON  Bx  TYPE  FIRST  STAGE  COMPRESSORS 

1.  These  operational  tests  were  done  on  two  occasions  altoge¬ 
ther,  once  in  1970  and  once  in  the  period  1975-76.  There  were  more 
than  20  individual  tests.  The  objective  of  the  tests  in  1970  was  to 
do  research  on  the  causes  of  the  blade  malfunctions  in  order  to 


select  a  design  for  the  new  1-02  type  blade  fop  use  with  the  fifth 
and  sixth  batches  of  engines.  The  experiments  done  in  1975  had  to 
do  with  work  that  had  to  be  done  before  finalizing  the  form  of  the 
type  engine  (01  batch).  Some  supplementary  experiments  were 
also  done  in  order  to  try  to  find  the  effects  of  installing  blades 
with  varying  frequency  differences  on  the  long-term  erosion  of  the 
blades  and  the  prevention  of  the  occurrence  of  vibration  in  order 
to  raise  the  operational  reliability  of  the  blades.  Besides  this, 
tests  were  also  run  on  the  similar  type  2C  Soviet  blade  for  purposes 
of  comparison. 

2.  The  contents  of  the  tests 

We  made  measurements,  when  S  —  0.7,  0.8,  0.85,  0.90,  0.95,  i.o  , 

of  the  overall  performance  of  the  compressor  and  its  base  character¬ 
istics  of  performance  when  n  *  1.0  as  well  as  of  the  surge  margins. 
In  order  to  carry  out  observations  of  the  blades  involved,  strain 
panels  were  attached  to  the  areas  in  which  malfunctions  often 
occurred  (see  Figure  2)  and  measurements  were  made  of  the  frequency 
of  the  movements  of  the  blades  and  of  their  changes  in  applied 
stress . 

3.  On  the  basis  of  the  fact  that,  in  the  factory,  there  were 
four  occurrences  of  cracking  and  one  occurrence  of  actual  breaking 
all  of  which  occurred  during  heat  test  configurations  and  outside 
of  the  factory,  there  was  one  case  of  bending  and  five  cases  of 
cracking  which  occurred  in  high  altitude  situations  in  which  large 

M  numbers  were  involved,  it  is  possible  to  estimate  that  its  reduced 
rotation  speed  is  n  *  0.825 — 0.86.  The  high  altitude  and  tempera¬ 
ture  reduced  speed  of  rotation  is  as  shown  in  the  table  below, 
that  is 


Prom  the  table  above,  one  can  know  that  n  =  Q.85  is  the  approx¬ 
imate  critical  rotational  speed  of  malfunction  and  that  this  speed 
requires  careful  watching  as  well  as  the  careful  selection  of  the 
stress  values  involved. 

4.  We  made  measurements  of  the  first  stage  performance  of 
the  compressors  being  tested  and  we  made  use  of  the  rotor  of  the 
same  test  bed  to  measure  n  and  G  for  the  mechanism  as  a  whole  and 
these  results  closely  followed  the  common  operational  curves  for 
the  compressors  and  the  turbines  involved. 

5.  New  topics 

On  the  basis  of  the  experience  of  our  factory  with  the  elimina¬ 
tion  of  malfunctions  as  that  experience  was  acquired  before  19.7CL, 
the  first  step  was  to  recognize  the  fact  that  blade  malfunctions 
were  generally  caused  by  vibration  fatigue  and  that  cracking  of  the 
blades  in  these  engines  took  place  when  M  *  0.85  and  the  engine  was 
in  a  heat  test  configuration.  Moreover,  such  an  appearance  of 
cracking  was  always  present  when  the  engine  involved  was  experienc¬ 
ing  abnormal  sonic  conditions.  When  the  period  of  time  Involved  in 
this  process  was  short,  it  was  possible  to  see,  in  the  blades  in¬ 
volved,  violent  agitated  movements.  On  the  basis  of  other  research 
materials  published  in  China  and  abroad,  we  came  to  realize  that  this 
type  of  phenomenon  was  stall  vibration,  that  is  to  say  that  in  cases 
where  the  phase  frequency  which  gives  rise  to  the  violent  mvoements 
is  the  same  as  the  phase  frequency  which  is  fixed  by  the  nature  of 
the  blades  involved,  one  gets  the  appearance  of  large  amplitude  vibra¬ 
tions  and  this  leads  rapidly  to  the  cracking  of  the  blades.  Also, 
due  to  the  fact  that  the  violent  vibrating  movements  of  the  blades 
causes  agitation  in  the  air,  one  gets  the  appearance  of  an  accompany¬ 
ing  whistling  sound. 

When  the  test  operations  were  begun,  three  of  the  rotors  were 
not  equipped  with  strain  plates.  In  the  process  of  running  tests  on 
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Figure  2.  Strain  plate  positions 

1 — top  section;  2 — 1  narrow;  3 — 2  thick;  4 — 3  narrow;  5 — 4  narrow; 

6 — 5  thick 

the  blades  from  the  fifth  batch  of  engines,  when  n  =  0.85,  one  got 
the  occurrence  of  an  unusual  whistling  noise  similar  to  the  whistle 
of  a  train  in  the  distance.  Later,  use  was  made  of  the  strain 
plates  in  order  to  make  observations.  Stresses  were  found  to  be 
as  high  as  450 0  kg/cm^.  Because  of  a  shortage  of  measurement  equip¬ 
ment  with  the  needed  capabilities,  we  were  capable  of  recognizing 
the  existence  of  vibration  states  of  the  blades  in  the  intake;  how¬ 
ever,  there  was  no  way  to  distinguish  the  different  natures  of 

these  vibrations.  For  the  results  of  these  tests,  see  Table  3  and 

Table  4  as  well  as  Figure  3- 

3.  AN  EXPLANATION  OF  THE  OVERALL  RESULTS  OF  THE  TESTS 

1.  For  rotors  in  the  fifth  batch,  except  for  the  case  in  which 

n  *  0.7  from  Q. 75-1.0,  one  sees  the  first  appearance  of  vibration 
points.  The  vibration  threshold  is  very  close  to  the  operational 
points  in  thi3  case.  This  is  particularly  true  for  the  case  in  which 
n  *  0.8-0.85  which  can  be  seen  in  Table  3  and  Table  4. 

2.  The  thick  and  thin  blades  in  the  fifth  batch  had  no  differ¬ 

ences  between  them  as  far  as  their  compositions  were  concerned.  The 
vibration  thresholds  for  the  former  were  shifted  somewhat  to  the 
left  as  compared  to  the  vibration  threshold  of  the  latter.  However, 
the  difference  was  not  obvious. 
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Figure  3.  Performance  of  B1  type  first  stage  compressors 

1 — sixth  batch  vibration;  2 — surge;  3 — sixth  batch;  A — surge;  5 — 
fifth  batch;  6 — vibration;  7 — sixth  batch;  8 — common  operational 
curve;  9 — fifth  batch;  10 — common  operational  curve;  11 — kg/sec; 
12— surge;  13— fifth  batch  surge;  14— fifth  batch  vibration; 

15 — kg/sec 


TABLE  4 


1-02 

Fifth  batch 

Rijcth  batch 

1C 

a 

G, 

n 

c. 

m 

G. 

C. 

Surge 

1.06 

49.09 

1.548 

1  .OF 

56.64 

l  .OF 

1  .OF 

50.5 

1.52 

vibra- 

0.9S 

44.26 

1.446 

0.95F 

53.76 

1.494 

0.95F 

1  .OS 

49.7 

1.52 

tion 

0.85S 

41.33 

1.393 

0.90F 

50.75 

1.447 

0.9F 

0.90S 

46.  ILF 

1.438 

thre- 

O.SS 

38. 19 

1.339 

0.85F 

45.41 

1.393 

0.9S 

0.90S 

45.5LGF 

1.438 

shold 

0.7S 

32.39 

1.25 

0.80F 

41.64 

1.341 

0.85S 

1.40 

0.80S 

39.4LF 

1.33 

0.70S 

32.51 

1.256 

0.80S 

40.42 

1.346 

0.80S 

39LGF 

1.33 

0.70S 

33.6 

1.239 

0.7S 

33.8LF 

34  LGF 

1.246 

mm-F  mm-s 


']'**  LF  t 


LGF 


1 — vibration;  2 — surge;  3 — small  frequency  difference;  4 — large 
frequency  difference 


3.  Concerning  the  surge  threshold  of  the  sixth  batch,  when 
n  *  0. 9-1.0,  first  there  is  vibration  and  then  surge.  The  amount 
of  flow  Jumps.  At  this  time,  AG  *  1.27  kg/sec.  At  other  speeds  of 
rotation,  there  is  surge  but  no  vibration.  The  boundary  is  located 
between  1-02  and  the  fifth  batch. 

4.  Concerning  the  boundaries  which  pertain  to  the  1-02  blades, 
there  is  a  distinction  that  must  be  made.  If  one  is  considering  the 
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case  In  which  there  is  a  large  frequency  difference  distribution, 
then  there  is  no  vibration  boundary.  In  the  case  in  which  there 
is  a  small  frequency  distribution,  that  is,  when  n  =  0. 9-1.0,  vibra¬ 
tion  does  occur.  In  the  tests  that  were  done  on  1-02  blades  in 
1970,  there  was  no  vibration. 

5.  Concerning  2C  blades,  they  are  basically  the  same  as  the 
1-02  blades;  however,  they  do  exhibit  vibration.  Prom  the  tests 
which  were  described  above,  it  can  be  seen  that  the  values  which  are 
selected  for  the  blade  density,  aspect  ratio,  chord  length  and  C^^ 
as  well  as  the  selection  of  varying  angles  of  attack  and  lag  angle 
do  not  Just  influence  the  amounts  of  flow  and  the  efficiency,  but 
they  also  influence  the  vibration  threshold  of  the  blades  involved. 
The  application  of  varying  frequency  differences  has  obvious  effects 
on  th  continuation  or  elimination  of  vibration. 

CONCLUSIONS 

1.  Thin  blades  from  the  fifth  batch  can  easily  exhibit  Intake 
vibration  when  n  =  O.8-O.85.  The  thick  blades  from  the  fifth  batch, 
when  n  *  0.8-0.85,  have  a  vibration  range  which  Is  somewhat  improved 
over  that  of  the  thin  blades.  2C  blades  exhibit  vibration  when 

n  *  0.95-1.0,  and  the  stress  on  them  is  2960  kg /cm  . 

As  far  as  a  comparison  of  the  three  types  of  moving  blades  is 
concerned,  the  1-02  blade  Is  capable  of  the  safest  operation. 

2.  Concerning  the  use  of  equipment  having  large  frequency  diff 
erence  distributions,  in  order  to  cause  an  improvement  in  the  vibra¬ 
tion  of  the  blades  being  considered,  the  use  of  such  a  large  fre¬ 
quency  difference  causes  the  boundary  line  involved  to  move  to  the 
left.  The  2C  and  1-02  blades,  because  of  the  small  frequency  dis¬ 
tributions  Installed  with  them,  still  do  exhibit  vibraiton.  However 
if  they  are  employed  with  large  frequency  differences,  then  it  is 
possible  to  eliminate  vibration.  The  author  knows  this  method  to  be 
effective. 
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3.  Looking  from  the  point  of  view  of  performance,  in  terms 
of  the  amount  of  flow  and  the  pressure  ratios  involved,  the  2 C 
blades  are  the  larger  and  the  1-02  blades  are  in  the  middle.  In 
terms  of  efficiency,  the  1-02  blades  are  the  highest,  the  2C  blades 
are  in  the  middle  and  the  fifth  batch  blades  are  slightly  lower 
than  the  2C  blades. 

4.  Experimental  error  as  reflected  in  a  comparison  between 
the  1970  and  1975  data. 

(1)  Error  in  amounts  of  flow  n  *  1.0  (when  the  throttle  is 

all  the  way  open,  for 
example ) 

1-02  AG  =  0.5  kg/sec 
fifth  batch  AG  =*  0.24  kg/sec 

2C  AG  =  0.1-0.33  kg/sec 

(2)  Concerning  the  pressure  ratio  error,  when  n  =  1.0  is 
taken  as  the  example,  the  1975  data  is  slightly  lower  than  the  data 
for  1970. 

(3)  As  far  as  errors  in  efficiency  are  concerned,  the  values 
of  nk  for  the  three  types  of  rotors  in  the  1975  tests  are  all  lower 
than  the  corresponding  values  from  the  1970  tests.  This  is  primar¬ 
ily  determined  by  measurement  errors  in  Gg,  and  N  (power). 

Considering  the  errors  described  above  as  a  whole,  the  tests 
which  have  been  spoken  of  in  this  article  had  to  do  with  the  deter¬ 
mination  of  the  displacements  and  direction  of  movement  of  the 
surge  and  vibration  boundaries  under  differing  conditions;  therefore, 
these  errors  can  have  no  substantial  effects  on  the  findings. 
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B,  TYPE  TRANSONIC  AXIAL-FLOWCOMPRESSOR 

ist  Stage  rotor  blade  design,  test  and 

APPLICATION 

Zhou  Wei-guo 

(Sknyyang  Aeromgiot  Company)  • 

Abstract 

Thi*  paper  describes  the  design  modification  of  the  1st  stage  rotor  blade  to  relieve 
crack  development  and  presents  test  result  of  the  modified  blades. 
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